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ABSTRACT 

The mid-infrared ratio [Ne m] i5.6/im/[Ne ii] n.%p_m is a strong diagnostic of the ionization state of 
emission hne objects, due to its use of only strong neon emission Hnes only weakly affected by 
extinction. However this ratio is not available to ground-based telescopes as only a few spec- 
troscopic windows are available in the MIR. To deal with this problem we aimed to verify if 
there exists a conversion law between ground-accessible, strong MIR line ratio [S iv]/[Ne ii] 
and the diagnostic [Ne iii]/[Ne ii] ratio that can serve as a reference for future ground-based 
observations. We collated the [Siv]io.5^m, [Ne ii]i2.8^m, [Neiii]i5.6^m and [S iii]i8.7^m emis- 
sion line fluxes from a wide range of sources in the rich Spitzer and ISO archives, and 
compared the [Neiii]/[Neii], [S iv]/[S iii], and [Siv]/[Neii] ratios. We find a strong correla- 
tion between the [Siv]/[Neii] and [Ne iii]/[Ne ii] ratio, with a linear fit of log([Neiii]/[Neii]) 
= 0.8 1 log([S iv]/[Ne n]) -1-0.36, accurate to a factor of ~2 over four orders of magnitude in the 
line ratios. This demonstrates clearly the ability of ground-based infrared spectrographs to do 
ionization studies of nebulae. 
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1 INTRODUCTION 

The mid-infrared spectrum (MIR, ~3-40//m) hosts a range of im- 
, portant diagnostic fine-structure emission lines. Due to the much 
' lower dust opacity at these long wavelengths, the emission lines 
' are significantly less affected by extinction than ultraviolet or opti- 
cal lines (see e.g. Draine 2003). Thus these lines are able to probe 
ionized regions that lie behind dense obscuring clouds or even deep 
within galaxy centers, giving insights into the regions where mas- 
sive star formation events actually occur, or heavily obscured active 
galactic nuclei (AGN) can exist. 

One of the more useful facets of this wavelength regime is that 
it includes emission lines from several ions of the same species, 
in particular strong lines from neon and sulphur. These lines can 
be used as strong diagnostics of the ionization state of the emit- 
ting gas as, being from the same species, there is no direct de- 
pendence on the gas phase abundances, and the different ioniza- 
tion potentials allow the observer to trace the hardn ess of the 
ionizing spectrum (see e.g. iDopita & Sutherland I2OO3I) . Lines ra- 
tios such as [Nev]i4.3^n,/[Nen]i2.8;jm, [Neiii]i5.6;,m/[Neii]i2.8,/m, and 
[S iv]io.5^ni/[S iii]i8 7^ni h ave been used to diagnose the presenc e of 
deeply buried AGN (e.g. lSturm et alj2002l : [Annus et alj2006l) and 
as indicators of the ages of starbursts (e.g. Achtermann & Lacvl 
ll995l : lThomlev et alj2()00l) . 

The major drawback of the MIR emission lines as diagnos- 
tics is their accessibility, with the full range of lines only avail- 
able to spectrographs aboard sp ace-borne telescopes lik e the short 
wavelength spectrograph (SWS. Ide Graauw et aljl996h aboard the 
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Infrared Space Observato ry (ISO, Kessler et lj] |l996h and the in- 
frared spectrogr aph (IRS. iHouck et alj|2004h aboard the Spitzer 
Space Telescope Jwemer et al. 2004) . While MIR instruments on 
ground based telescopes can provide higher spatial and spectral 
resolution, only a sparse number of MIR bands are available due 
to the opaqueness of the Earth's atmosphere in the mid-infrared 
(M-band:~ 4.5 - 5pm, N-band:~ 8.5 - I2>pm, and Q-band:~ 17 
- 23/jm). As a result, the [Neiii]i5 5,,n, line is inaccessible from the 
ground and therefore recent high-resolution studies h ave resorted to 
the ra tio [S iv]io.5;jm/[Neii]i2.8;jm as a diagnostic (e.g. lSniiders et al.l 
|2007|) . in a similar manner to earlier aircraft-borne IR telescope 
studies (e.g. I Achtermann & Lacvl[l995 l). However the question re- 
mains whether the use of this ratio as a proxy for the ionization sen- 
sitive [Neiii]i5 6,,n,/[Neii]i2.8,ira fatio is justified. While photoioniza- 
tion theory clearly indicates that there will be a correlation between 
these ratios, the dependence of the line emission on often unknown 
physical parameters make the accuracy of this correlation uncer- 
tain. As an example, the [S iv]/[Ne 11] ratio is directly proportional 
to the sulfur to neon abundance ratio, while the [Ne iii]/[Ne 11] will 
be almost unaffected by any changes in the relative abundances. 
Similarly, with an excitation potential of 34.79 eV, [S iii] is sensi- 
tive to lower energy photons than [Neiii], which has an excitation 
potential of 40.96eV and will not exist in photoionized gas when 
there are no photons above this limit. 

In this letter we use archival ISO and Spitzer observations of a 
range of astrophysical objects to demonstrate that such a correlation 
does exist and derive a simple conversion law for these ratios. We 
mention possible systematics with the determination of these ratios 
and discuss the origin of the observed correlation. As a final note 
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Table 1. Publications from which we sourced our MIR emission lines dif- 
ferentiated the object classes, along with the telescope used and the mean 
fractional eiTor in the [Neii]i2.8fira line- 



Objects 


o"Ncll/[Neii] 


Publication 


• 

Galactic H ii regions 


n in 


Lriveon et al. ( 2UU2 ) 


vjlclill n. 11 IC^A'-'AiS 


15 


T phniitpillpr pt al (')(M^^^^^ 


Extragalactic 




H II regions & nuclei 


0.05 


Dale et al. (20Q6)' 


MlOl Hii regions 


0.39 


(jordon et al. (2008)' 


M83 H II regions 


0.01 


Rubin et al. (200T)' 


M33 Hii regions 


0.01 


Rubin et al. (2008)' 


M82 Pointings 


0.15 


Beirao et al. (2008)' 


Arp244 Pointings 


0.04 


Brandl et al. (2008)' 


Starburst galaxies 


0.09 


Brandl et al. (2006)*'' 


Starburst galaxies 


0.30 


Verma et al. (2003)^ 


BCDs 


0.04 


Wu et al. (2008)' 


BCDs & Starbursts 


0.09 


Enselbracht et al. (2008)' 


ULlRGs 


0.10 


FaiTah et al. (2007)' 


AGN 


0.42 


Sturm et al. (2002)^ 


AGN 


0.14 


Deo et al. (200T)' 


LMC & SMC PNe 


0.11 


Bemard-Salas et al. (2008)' 



1 : Spitzer 

2 : ISO-SWS * : Private Communication 



we recommend the application of these results in future (ground- 
based) observations. 



2 OBSERVATIONAL DATA 

The strong, diagnostic, MIR emission lines [S iv]io.5;,m, [Neii]i2.8//m, 
[Neiii]i5 5,,n,, and [Sm]ig7,,n, have been detected and characterised 
in a broad range of astrophysical objects, from nearby galactic H ii 
regions to ultraluminous infrared galaxies (ULIRGs). To determine 
the connection between the [Siv]/[Neii] and [Nein]/[Nen] ratios 
we have extracted from the recent literature the fluxes of these four 
emission lines. While not exhaustive, the sample of 355 emission 
line objects is representative of the range of sources from which 
these lines arise, covering a range of physical parameters, such as 
source morphology and geometry, nature of the ionizing sources 
and metal and dust abundances. 

A list of the difi'erent source types in our sample and the 
relevant papers from which we extracted these is given in Table 
[T] Included in the sample are sources ionized by massive stars, 
both small scale (Hn regions and extragalactic Hii regions) and 
large scale (blue compact dwarfs (BCDs), starburst galaxies (SB), 
and Ultra Luminous IR galaxies (ULIRGs)), planetary nebulae 
(PNe), which are ionized by white dwarfs, and active galactic nu- 
clei (AGN). 

The table also includes the average percentage uncertainty in 
the [Nen] line as given in the reference. This uncertainty varies 
greatly between the various references, as the difi'erent authors have 
used different methods to account for both the observational and 
systematic uncertainties for their sources, which vary in flux and 
exposure times. While these uncertainties may bias our fit, the close 
correlation observed and the varying sizes of the samples (both dis- 
cussed in the next section) will tend to minimize these effects. The 
given line uncertainties in the references were carried over when 
calculating their corresponding ratios. 

The majority of objects within our sample were observed with 



the Spitzer IRS, with the remainder coming from ISO SWS ob- 
servations (as noted in Table [T]!. In the cases where objects were 
observed twice (for example with both Spitzer and ISO), we have 
included both observations, to indicate the systematic uncertainty 
in the sample. For all of our sample, all four MIR lines fall within 
Spitzer s IRS short-high (SH) module (9.9-19.6yum) and ISO's 
SWS (2.38 to 45.2/jm), meaning that aperture effects are limited 
to that due to the difi'erent wavelengths of the lines in the ratios (a 
factor of 1.2 for [Siv] to [Neii] and [Neii] to [Neiii]). 

Within our total sample there are approximately 8 objects ob- 
served by both telescopes, and a similar number of objects observed 
by Spitzer whose lines have been extracted by two difi'erent authors 
using different methods. For objects by both telescopes (e.g. NGC 
4151), the difference in line ratios is of the order of 0.3 dex, most 
likely due to the different telescope apertures and pointings. Some 
objects were observed once by ISO, but have multiple pointings by 
5p//zer-IRS (e.g. M82), and for these objects we consider the ISO 
observation an integrated fiux and the IRS observations to be indi- 
vidual Hii regions. For the few IRS observed objects investigated 
by different authors, the differences in line ratios arise purely due to 
the different line fiux extraction techniques (e.g. IIZw40) and are of 
the order of 0.2 dex, and within the given uncertainties. However 
the small number of objects observed multiply (~15, distributed 
across both AGN and Starburst/BCD) compared to the total sample 
size (355), mean that a bias of the final fit towards these types of 
objects will not occur 

In extracting our sample we excluded all objects in which any 
of the four lines were undetected, or for which only upper limits 
were provided. This means that we do have a bias in the sample 
against weak line objects, but, with objects covering a broad range 
of distances and four orders of magnitude in line ratio, this bias is 
likely to be weak. This bias is present, however, in at least one of 
the groups (ULIRGs), and is discussed further in the next section. 

The total sample represents a large variety in ionising condi- 
tions (-2.1 <log([Neiii]/[Neii])< 1.8), metal abundances (ranging 
from log(0/H) -I- 12 « 7.1 (SBS 0335-052W) to « 8.8 (IC342)) 
and the object particulars such as infrared luminosities (PNe to 
ULIRGs), galaxy types (blue compact dwarf galaxies to massive 
ULIRGs), and star formation rates (BCDs and Hii regions of less 
than a few yr"' to ULIRGs ~ lOO's yr"'). 



3 RESULTS & DISCUSSION 

The final sample of 355 emission-line objects is presented in figure 
[T]in the form of two emission line diagnostic diagrams: [S rv]/[S m] 
vs [Nein]/[Nen] (figure [T)i) and [Siv]/[Neii] vs [Ne iii]/[Ne ii] (fig- 
ure [TJj). The separate classes of object within the sample are dif- 
ferentiated both by colour and shape, as labelled in the key in the 
upper left hand comer of the diagrams. In the lower right corner of 
the diagrams we show the average observational uncertainty across 
the sample, as given within the references. As noted previously, the 
uncertainty given in the references varies significantly between the 
different classes, ranging from -1% for M33 & M83 to -40% for 
the galactic Hii regions. In all classes however, the uncertainty on 
the [Ne iii]/[Ne ii] ratio is lower than that of the other two ratios 
considered here, partly due to the strength of the neon lines (higher 
signal-to-noise), but mostly due to the location of the [Sivjio.s^m 
in the silicate absorption feature (discussed further below), which 
both decreases the strength of this line and adds an additional un- 
certainty due to extinction. In any case, these uncertainties are ap- 
plied to the reduced x' fit discussed below. 
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Figure 1. Mid-IR line ratio diagrams revealing the correlation be- 
tween [Ne III] i5.6,,m/[Nen] 12.8pm and [S iv]io.5pm/[S iii]i8.7pm (a), and 
[Ne III] i5.6p,^/[Nen] 12.8pm and [S iv]io.5pm/[Neii] 12.8pm (b). Each diagram 
contains the observed ratios from several different emission-line objects, 
differentiated by colour and shape, as labelled by the key in the upper left 
hand comer, and discussed in section |2] The en'or bars in the lower right 
of the diagrams indicate the average uncertainty for the whole sample. The 
arrow indicates the effect of extinction on the ratios for an Ay = 30 mag- 
nitudes. The lower diagram also shows our best fit to the ratios (solid hne) 
along with the l-o" dispersion (dotted lines), as given in Table|2] 

In both diagrams we plot a vector that indicates the effects of 
an Ay = 30 magnitudes of extinction on the lin e ratios by a uni- 
form dust screen, assuming a lLi&Drainel ( l200lh extinction curve. 
Note that, depending on which curve is assumed, the length and 
direction of the extinction vector can change but that in all cases, 
the [S iv] 10.5pm line suffers the highest relative extinction (see Table 
2 and associated discussion in Farrah et al. 2007). 

The first diagram, figure [T^, reveals the standard space-based 
diagnostic for the average ionization state of the emission line ob- 
ject. This diagram is a strong diagnostic as it uses 4 strong lines, 
with each ratio involving two different ionization states of the same 
element, meaning that abundance differences are not an issue. Due 
to this ionization dependence and insensitivity to abundance, this 
diagram has been often used in the literature as a measure of the 
average hardness of the ionizing radia t ion fie ld, most recently by 
iGordon et al] feOOSh . In iGordon et all j2008l) , the correlation be- 



Table 2. Fit parameters for the full sample and individual classes, including 
the number of objects (N) within the sample. 



Object class 


N 


a 




o" r\x 


Full sample 


355 


0.81 


0.36 


0.25 


H n regions 


65 


0.74 


0.42 


0.23 


Extragalactic H n 


97 


0.82 


0.36 


0.15 


ULIRGs 


20 


0.49 


0.26 


0.33 


AGN 


20 


0.69 


0.29 


0.25 


PNe 


20 


0.74 


0.56 


0.30 


Starburst 


56 


0.65 


0.32 


0.28 


BCDs 


49 


0.86 


0.32 


0.30 



tween the [S iv]/[S iii] and [Neni]/[Nen], clearly visible in our fig- 
ure [T^, was quantified for both individual Hn regions in MIOI 
and for a combined MIOI and starburst/BCD galaxy sample from 
Engelbracht et al. (2008) (their equations 5 and 6 respectively). 
They find a steeper slope than appears in figure [T^ (1.53 versus 
1.29), but thi s may be due to th e smaller sample size, and is also 
biased by the lEngelbracht et"aD {2OO8) starburst sample. This sam- 
ple, marked as BCDs in our sample, can be seen as the offset empty 
circles at low [Ne iii]/[Ne n] in figure[T^. The visible offset from the 
rest of the sample may be intrinsic, but may also be from the un- 
derestimation of the contr ibution from the I 2.7/jm PAH feature in 
these weaker line objects jSmith et al.ll2007l) . 

There are several offset groups within figure [TJi, the most no- 
ticeable being the AGN. This group is clearly offset towards higher 
[S iv]/[S m] values relative to their [Ne in]/[Ne n] ratio. The reason 
for this offset is uncertain, but is most likely related to the harder, 
power-law ionizing spectrum found in AGN. This leads to higher 
average electron temperatures in the S^^ and Ne+ zones and differ- 
ent extents, relative to gas ionized by stars, and therefore a weaker 
[S in]/[Ne 11] ratio for AGN relative to stellar photoionization. This 
olTset may also be measurement related, with the weaker 12.7/jm 
PAH feature meaning that the [Neii] line is determined more ac- 
curately and with higher flux. The presence of ULIRGs and some 
starburst galaxies in this region as well is unsurprising as some of 
these sources contain known active nuclei, just as many of the AGN 
classified galaxies also have strong, active star formation. 

In figure [TJ) we show the ground-accessible ratio [S iv]/[Ne 11] 
versus [Neni]/[Ne 11], revealing a strong correlation between the ra- 
tios, justifying the use of the [S iv]/[Ne 11] ratio as a surrogate for the 
[Ne ni]/[Ne 11] ratio when unavailable. Compared with figure[T^, the 
data show surprisingly less scatter. This probably arises due to the 
use of a common denominator ([Ne n]), but the [S iii] line may also 
be introducing some scatter into figure[T^, possibly due to its longer 
wavelength (and hence lower resolution). 

To quantify this correlation we perform a linear least-squares 
fit to the full sample using the IDL routine FITEXY, giving the 
empirical relationship, 

log([NeIII]/[NeII]) = or log ([SIV]/[NeII]) (I) 
The result of this fit, along with the 1 cr dispersion, is given in Table 

m 

The first thing to note is that the log-space linear fit to the full 
sample appears to be accurate to 0.25 dex (that is, less than a fac- 
tor of ~2 in the linear ratio) over four orders of magnitude in the 
line ratios, demonstrating clearly the ability of ground-based IR 
spectrographs to do ionization studies of nebulae. The dispersion 
around this relation is greater than the average uncertainty, indi- 
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eating either inaccurate uncertainties, intrinsic variations due to the 
variety of sources or, most likely, a combination of both. 

In addition to the full fit we also present fits to our various ob- 
ject classes, also presented in Table|2l Overall, the relations match 
reasonably well with the full sample fit, with the two major out- 
liers either side being the planetary nebulae (high offset in p) and 
the ULIRGs (shallow slope in a). These individual fits can be used 
when the specific class of objects being observed are known, but 
we suggest that the full fit should be used in most circumstances, 
as it is a more robust fit and includes all object classes. 

In figure [TJ) the PNe are visibly offset from the rest of 
the sample, both to high [Siv]/[Neii] and higher than average 
[Ne iii]/[Ne ii]. This offset is most likely due to the hotter and harder 
white dwarf ionizing source in PNe, although this then poses the 
question why the AGN, which also have a harder ionizing spec- 
trum than that found in H ii regions and galaxies, is not also more 
offset. 

The shallow slope of the ULIRG sample demonstrates one of 
the issues with the [S iv]/[Ne ii] ratio: the locati on of the [S iv]io 5um 
in the 9.7/[jm silicate absorption feature. Of the lFarrah et al.l ^2007 ) 
sample of 53 ULIRGs over half have no detections or upper limits 
for the [S iv] line, while a large fraction of those which do have 
detections are uncertain, with detections of less than 3 sigma. As 
noted in the paper, only in ULIRGs with low silicate depths (5 sii 5 
2.1) is the [Siv] feature detected, even when the [Neiii]i5 6,,n, line 
is present. Thus there is likely both a bias in our ULIRG sample 
and an offset, such that low [Siv]/[Neii] ULIRGs are also those 
with high extinction and therefore offset from the relation, while 
the highest extinction sources are not included at all.. 

Another spectral feature which may affect the line ratio rela- 
tion is the 12.7^m PAH feature seen clearly in many galaxy spectra 
( ISmithetalJl2007l) . This broad feature lies below the [Neii] line, 
and makes the determination and subtraction of the underlying con- 
tinuum difficult, and thus increases the uncertainty in the flux of this 
line. As mentioned befo re, this may explain the offset seen in the 
lEngelbracht et all ( l2008h sample. 

However, even with these outliers the correlation is surpris- 
ingly tight. Simple photoionization theory expects such a correla- 
tion, with higher ionization leading to higher values for both ratios 
jPop ita & Sutherland 2003), but variations in the ionization param- 
etexJ and ionizing spectrum cause different changes in the two ra- 
tios. Such variatio ns have been explored in th e case of starbursts 
and Hn regions in iThomlev et al] bOOOt) and iRubin et"ai] bOOSi) 
and, given these results, the tight fit indicates a close correlation 
between ionization parameter and spectral hardness in these ob- 
jects. 

One further issue is abundance variations. The [Ne iii]/[Ne ii] 
ratio is sensitive to abundance variations only indirectly, through 
the resulting temperature effects. The [S iv]/[Ne ii] ratio however, 
can be directly affected by variations between the sulphur and neon 
abundances. As both are primary elements, significant total abun- 
dance variations are not expected. However, depletion of sulphur 
onto dust may cause variations in the gas phase and hence in the 
emission lines. While the total depletion of sulphur onto dust is un- 
certain (though see .Leboiiteiller et al. 2008) , it does not seem to be 
a gre at effect on the Ne/S ratio, as shown by iRubin et alj j200l 
l2008h . and therefore on the [S iv]/[Ne ii] ratio. 

Thus, in summary, to determine whether it is reasonable to 



The ionization parameter is a measure of the ionizing photon density over 
the gas density; U = gtot/"HC, where 2tot = Jis ^■'//'v dv 



use the ground-accessible ratio [S iv] io.5,im/[Ne ii] i2.8;jm as a replace- 
ment for the diagnostic [Ne iii]/[Ne ii] ratio, we have collated exist- 
ing ISO and Spitzer spectral observations of a wide range of as- 
trophysical objects from the literature. We find a good correlation 
between these ratios with a linear fit giving the relation, 

log([NeIII]/[NeII]) = 0.81 log ([SIV]/[NeII]) -I- 0.36, 

with a 1 (T dispersion of 0.25, corresponding to an uncertainty in the 
estimated line ratio of ~ 70%. We propose that for future ground 
based observations this relationship be used to determine the ion- 
ization state of the observed objects, though note that caution must 
be applied when looking at heavily extinguished (Ay S 20) objects 
where the [S ivJio 5,,„i may be affected by the 9.7/jm silicate absorp- 
tion feature. 
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